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We have developed a formalism suitable for calculation 
of the output spectrum of a detector continuously measuring 
quantum coherent oscillations in a solid-state qubit, start- 
ing from microscopic Bloch equations. The results coincide 
with that obtained using Bayesian and master equation ap- 
proaches. The previous results are generalized to the cases of 
arbitrary detector response and finite detector temperature. 



I. INTRODUCTION 

Quantum coherent (Rabi) oscillations in a two-level 
system represent a simple ancLfundamental example of a 
nontrivial quantum behavior.EI Recently increased inter- 
est in this subject is obviously related to the use of two- 
level systems (qubits) as, building blocks of a prospec- 
tive quantum computer.0 The emphasis has naturally 
shifted from traditional observations of Rabi oscillations 
in ensembles of two-level atoms to the studies of sin- 
gle qubits. Concentrating in this paper on solid-state 
qubits, let us mention recent demonstrations of single- 
qubit quantum coherent oscillations in both timeo and 
frequencjoB domains. (Rabi oscillations in individual 
quantum dots have been also demonstrated by traditional 
optical means;El however, we will discuss only solid-state 
qubits with electronic readout.) _ _ 

Even though experiments of Refs.EfEl have been done 
with single qubits, their results are essentially ensemble- 
averaged, and as a consequence, the problem of a quan- 
tum state collapse due to measurement is not very im- 
portant. Another possible experimental setup (for which 
the collapse is of the major importance) is a continu- 
ous monitoring of single-qubit quantum oscillations (Sg- 
|l|). Such setup is similar to that used in experimentsQ'H, 
with the difference that the detector output signal /(t) 
was actually not studied experimentally. 

The basic questions in prospective experiments of the 
type shown in Fig. |^ are the following: 1) what is the 
effect of continuous quantum measurement on the qubit 
evolution, 2) how the detector output I{t) looks like, and 
3) what is the relation between detector output and qubit 
evolution? Some answers to these questions have been 
obtained recently (see, e.g., Ref.a and references therein); 
however, the subject is still active and quite controversial. 
(The main reason of controversy is the necessity to go 
beyond the Schrodinger equation to describe continuous 
collapse due to measurement.) 

In this paper we will consider a relatively simple ques- 
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FIG. 1. Schematic of a single solid-state qubit continuously 
measured by a detector. 



tion: what is the spectral density Si{uj) of the detector 
output I{t) and how high is the spectral peak correspond- 
ing to quantum oscillations of the qubit state? (Notice 
that we assume detector output I{t) to be a classical 
magnitude, so Si{uj) does not depend on a particular 
method of further signal processing.) This question has 
beenr8ud|djessed already in a number of papers (see, e.g., 
Refs.tHila) using various techniques. 

Spectral densityrj^Sf (w) has, been calculated using 
Bayesian formalismQ'Ej in RefJl3 for the case of a weakly 
responding (hnear) detector. In particular, it has been 
shown that the spectral peak Si (ft) at the frequency 
of quantum oscillations cannot be higher than where 
5*0 is the noise pedestal due to intrinsic detector noise. 
It has been also shown that the results of the BayesiaH. 
formalism for Si{u)) exactly coincide with the resultaij'Ell 
obtained using the standard master equation formalism. 
In spite of the same results, the interpretations are quite 
different since the Bayesian formalism describes individ- 
ual monitoring of quantum evolution in time and treats 
I{t) as a classical measurement result, while the master 
equation formalism can describe only ensemble-averaged 
magnitudes and in some sense should treat I(t) as a quan- 
tum operator. (The results of two formalisms coincide 
because Sj (uj) is essentially an average quantity. In more 
general |aituations, for example, for a quantum feedback 
analysiajEiJ the results are not the same since the master 
equation formalism fails.) ■— . 

The results for Si{uj) have been confirmed in Ref.Ej 
using a somewhat different approach based, on the gen- 
eral theory of linear detectors. In Ref.Ej the results 
have been pCaniirmed using the approach of quantum 
trajectorieScj adopted from quantum optics (this ap- 
proach is siiailar to the Bayesian formalism) . It has been 
also shownEj that the same formulas for Si{uj) remain 
valid even when the condition of weakly responding de- 
tector is not satisfied (the detector response A/ to the 
change of the qubit state is comparable to the average 
current Iq). j-, 

A different result for Si{uj) has been obtained in Ref.EJ 
(the studied system is slightly different from the type 
shown in Fig. ^ however, it is essentially similar). In 
particular, the calculated ratio [<S'/(ri) — So]/So can be 
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arbitrary large (not limited Jay 4). In our opinion, this 
is because the result of Ref.tJ includes the contribution 
from zero-point oscillations which are not measurable in 
a straightforward way. In principle, zero-point oscilla- 
tions can be sensed in an experiment; however, it would 
necessarily require quantum (non-classical) interaction 
between the detector [which outputs I{t)] and the next 
stage. In other words, it would require measurement of 
Si{(jj) without ever measuring I{t). A simple example 
of such setup is the absorption/emission of photons at 
resonant frequency fi. A more sophisticated example of 
such measurement of qubit oscillations has been consid- 
ered in Ref.li3. The idea is to use rotating measurement 
basis, in which there are essentially no oscillations, but 
rather jumps between two stationary states due to exter- 
nal noise. Formally shifting the zero-frequency spectrum 
of such continuous measurement to the frequency Q, one 
can obtain arbitrary high spectral peak. (There is no re- 
strictie]*-on the ratio Si{lo — 0)/So in a strong dephasing 
case.t2rilj) Such setups, however, are not the subject of 
the present paper since we limit ourselves to the straight- 
forward case of Fig. I with I{t) being usual classical sig- 
nal which can be amplified further by any good amplifier. 
Also, we do not consider here the quantum feedback se- 
tups, which can provide arbitrary high spectral peak of 
the detector current at the oscillation frequency ri.E5l 

The main result [5/(51) — So]/So < 4 for a simple setim 
seems to contradict the experimental results of Ref.E^ 
which claim the measurement of Si{uj) from a single spin 
precession in an STM-based setup (a significantly dif- 
ferent, but still an spalogous experiment). In a simi- 
lar recent experimentcj the maximum observed peak-to- 
pedestal ratio for a measurement of a single spin preces- 
sion was a little less than 4x3 We cannot explain the dis.- 
agreement between the theory and experiment of Ref.Ei 
however, we note that in the recent theoretical paperES 
which considers a somewhat similar setup, the possibility 
of a relatively high spectral peak has not been confirmed. 
The spectral peak of Si{uj) due to qubitj-ascillations has 
been also considered theoretically in RefO; however, the 
peak magnitude has not been calculated. 

Let us finally mention one more theoretical approach 
(developed by S. A. Gurvitz) to the analysis of detec- 
tor output I{t) based on Bloch equations which describe 
the evolution of the coupled "qubit -|-detector" density 
matrix.t3 The advantage of this approach is the straight- 
forward microscopic derivation of Bloch equations from 
the Schrodinger equation, while in the approaches men- 
tioned above (Bayesian, master equation, quantum tra- 
jectory) the collapse ansatz is either explicitly or implic- 
itly used. The Bloch equation approach has been usedJx) 
calculate some statistical characteristics of I{t) in Ref.EEl; 
however, the spectral density Si{ui) has been obtained 
only at the frequencies lower than the ensemble dephas- 
ing rate and so the spectral peak due, to quantum oscil- 
lations has been out of scope of Ref£3. 

In this paper we show how the Bloch equations can 
be used to calculate the detector output spectral density 
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FIG. 2. Small transparency QPC (tunnel junction) as a 
qubit state detector. The barrier height depends on the elec- 
tron position in the double-dot system. 



Si{w) for a particular measurement setup (Fig. g). The 
case of an arbitrary qubit coupling with detector and fi- 
nite detector temperature is considered. We prove that 
the results for Si{lu) coincide with that obtained previ- 
ously (in a narrower validity range) by the master equa- 
tion and Bayesian approaches. The Bayesian results are 
generalized to the case of arbitrary response factor and 
finite detector temperature; it is shown that the equiva- 
lence of results of the three approaches still holds in this 
case. 



II. THE SYSTEM AND BLOCH EQUATIONS 

We consider the system (Fig. ^ irgrptUi<iedj-H^^ 
and studied extensively after that.QT^EIIoSES The 
qubit is represented by a single electron in a double quan- 
tum dot. The detector is a quantum point contact (QPC) 
whose barrier height depends on the electron position, so 
the current / through QPC measures the qubit state in 
the basis of localized states |1) and |2). We will limit 
ourselves by the case of small transparency QPC which 
is equivalent to a simple tunnel junction. The Hamilto- 
nian of the system, 

= Hqb +'HdET +'HlNT, (1) 

describes the qubit, detector, and their interaction: 

,t.„ , J 

'Hdet = ^ Eia\ai + ^ Erolar 

I r 

+ ^Ma|a^ + H.C., (3) 

l,r 

TiiNT = ^ At- (4^1 - 4^2) ajar + H.c. (4) 



(for simplicity M and AM are assumed to be real and 
energy-independent). The average detector currents cor- 
responding to the qubit states |1) and |2) are equal to 
h = 27r(M + AM/2fpipre^V/n and h = 2n{M - 
/S.M/2Y piPre^V/h, correspondingly iV is the voltage 



'X-QB = I {c\ci - 4c2) + H {c\c2 + 4ci), (2) 
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across the tunnel junction, e is the electron charge, and 
pi^r are the densities of states in the electrodes), while 
the corresponding detector noises have white spectrum 
and are given by the Schottky formula: 



5*1,2 = 2e/i,2. 



(5) 



Note that the detector voltage V is assumed to be large 
enough, so that the typical quantum noise freauency 
eV/fi is much higher than all relevant frequencies .EZl 

In the following we will distinguish the weakly- 
responding limit, I A/ 1 ^ /q, where A/ = Ii — I2 is the 
detector response and Iq = (/i + h)/'^, and the finite 
response case, |A/| ~ Iq. Notice that the word "cou- 
pling" is reserved for a different combination of param- 
eters: C = h{Mf/SoH [here Sa = {Si + S'2)/2], which 
affects the quality factor of quantum oscillations of the 
qubit. The frequency of unperturbed oscillations (with- 
out detector) is equal to = \^4H^ + e-^/h, where H is 
the qubit tunneling matrix element (assumed to be real) 
and e is the qubit energy asymmetry. |_. 

Our starting point is the Bloch equationsEEl describing 
the ensemble averaged evolution of the density matrix pf^ 
in which the subscripts = 1,2) label the qubit state 
while n is the number of electrons passed through the 
detector (only diagonal in n matrix elements are consid- 
ered because the nondiagonal elements decay very fas- 
For our system the Bloch equations are the following^ 

Pll— Pll ^ Pll ^ Pll 

e e e 



P22 — ' ~ P22 ^ ~ P22 ^ ^ ~ P22^ 

e e e 

+ 2|lmp5'2. 



(6) 



(7) 



_ It + II + It + I2 ,n , V ^1 ^2 ^„_i 



P12 



2e 



■Pi2 



■P12 



I\ I-) F H 

n+l I • ^ n I • / n n \ /n\ 

P12 +l^Pl2 +1— (Pll -P22)- (8) 



h 



Here 



1 ~ exp{~eV/T) 



I- = 1+ eM-eV/T) (9) 



are the partial currents in two directions {Ii — 1^ — I^) 
and T is the detector electron temperature. Notice that 
this tcrrmprature is different from what was considered in 
Rcfs£3 Ej. In those papers the effect of nonzero temper- 
ature of a passive environment coupled to the qubit was 
studied, so the important parameter was T/hfl. In this 
paper (similar to Ref.E3) we consider the effective detec- 
tor temperature and the important parameter is T/eV, 
while a finite temperature T always implies T 3> hfl. 
(The phonon temperature in the vicinity of the qubit may 



still be low, since T is only the electron temperature.) 
The density matrix pf^ obeys the natural normalization 
condition E„(pl\ + P22) = 1- n a. 

Notice that tracing Bloch equations (^)-(||) over the 
detector degree of freedom n, one obtains a conventional 
master equation for the qubit: 

Pll = -2 — Impi2, pii+P22 = l, (10) 
n 

s H / 

P12 = i ^ P12 + i-^ (pii - P22) - rpi2, (11) 

where pij =^L-P,j- and the ensemble decoherence rate 
r is equal tocStJ 



ir- 



it- 



2e 
2e 



2e 

coth(ey/2T) . 



(12) 



(The decoherence rate 7 for a single system without en- 
semble averaging is different - see Section V.) 



III. SPECTRAL DENSITY VIA MACDONALD'S 
FORMULA 

The Bloch equations couple the qubit evolution and 
the number n of electrons passed through the detector. 
So, to calculate the spectral density Sj{uj) of the detector 
current, we need to express Si{ll}) in terms of n. This can 
be easily done for the clas^cal random process I{t) using 
the MacDonald's formulaEj 



Si{uj) r sill{ujT)dT 

Jo d-T 



(13) 



where {Q^{t)) = I{t')dt' - (/)t) ) and averaging 

is over Jime t (MacDonald's formula have been also used 
in Ref.E3). In our case the average current is equal to Iq 
(see below) for a non-zero qubit tunneling H, so 



{Q'{T))=e'{n\r))^{Iory^ 



(14) 



where {11? {t)) is the average square of the number of elec- 
trons passed through the detector during time interval r. 

To calculate (^^(t)) we can use the Bloch equations 
and the obvious relation 



{n^{r))^yy[Pii{r)+pUr)]- 



(15) 



However, the situation is not too simple because the left 
hand side contains the averaging over time while the right 
hand side is essentially ensemble averaging which de- 
pends on the initial condition p^ (0). Quite naturally we 
should assume n = at r = 0, so that p" (0) = (5„oPij(0), 
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but the question about the choice of Py(0) remains un- 
clear because the qubit state actually oscillates in time 
(for a nonzero H, which case we always assume be- 
low). A natural choice is to use the stationary value: 
Pij^st = limt_,oo I]n[P"i(*) + Pij^^)]^ ^^'^ it is possible to 
prove that this choice is really correct in the following 

As we know from the Bayesian formalismjj we can 
monitor the oscillating evolution of the qubit density ma- 
trix Pij{t) in an individual realization of the experiment 
using the detector output I{t). This at least means that 
Pij (t) exists (even though it cannot be obtained using 
Bloch equations because they imply ensemble averaging) . 
So, the correct procedure of calculating (rt^ (r)) would be 
the following. The right hand side of Eq. (|lj) should 
be calculated for various initial values Pijir = 0) corre- 
sponding to values pij (t) in a sufficiently long realization 
of a process, and then the result should be averaged over 
the time t [i.e. weighted proportionally to the occurrence 
frequency of various pij]. Now it is very important that 
the Bloch equations (||)-(||) are linear in respect to the 
initial condition. This means that instead of averaging 
the result for ri^ir) over initial condition pijij = 0), we 
can use the initial condition which is itself the value av- 
eraged over time, i.e. stationary value Pij.st discussed 
above (of course, we implicitly use the process ergodic- 
ity). This ends the proof. 

Thus, to calculate Si{lo) we should solve the Bloch 
equations starting from the stationary initial condition 
P"j(0) = SnoPij,st, then calculate (^^(t)) using Eq. (|l5), 
and then use MacDonald's formula (O) to obtain Si{uj). 
Notice that the stationary state pij_^t can be easily ob- 
tained from the master equations (p^)-(pT|) and the con- 
dition Pij = 0, that gives (at H ^ 0) 



z{r) 



EKi(^) 



(20) 



Pll.st = P22,st = 1/2, Pl2,st = 0. 



(16) 



(The stationary state would be different if ilae-nubit had 
an extra coupling to a passive environment JiJO however, 
we do not consider such case.) 

One can use this method to calculate Si{uj) in a 
straightforward way (we have done it numerically); how- 
ever, it is better to use an analytical simplification cal- 
culating directly d{n'^ {t)) / dr = E„ ^^^I/Onl'^) + P22('^)]- 
Using Eqs. (^ and (R) and shifting summation over n in 
terms containing , one gets the equation 



d{n\T)) ^ lo 
e 

AI 



dr 



2(n(r)) 
A{t) + 



coth 



eV 
2T 



z(t) , eV 
coth — 
2 2T 



where 



^(r)^EnKi(T)-p^2M], 

n 

(n(r)) =EnKi(r)-fp^2M] 



(17) 

(18) 
(19) 



Notice that z(t) = since the evolution starts from the 
stationary state, pfj{0) = SnoPij,st, so the corresponding 
term in Eq. (|l^) vanishes. 

To calculate (n(r)), we again use Eqs. (|) and (§, shift 
summation over n, and obtain the equation 



d(n(r)) /dr ^ Iq + z{t)AI/2. 



(21) 



Since the last term vanishes because of z{t) = and since 
(n(0)) = 0, we obtain a simple result {n{T)) = Iqt. In 
particular, this means that the average detector current 
is equal to Iq (we have used this result above). 

One can see that the term 2/o(n(r))/e = 2/q^t 
from Eq. (|l^) exactly cancels the contribution from 
the derivative of the last term of Eq. ([T^). The term 
(/o/e)coth(ey/2r) from Eq^^(|l^ after being plugged 
into MacDonald's formula (|l^) gives the constant noise 
pedestal 2e/o coth(ey/2T) (as usual, we should use a 
smooth integral cutoff at high frequencies). In this way 
we get equations 

Si{uj) = So coth h 2weA/ / A{t) sin(cjT) dr 

27" Jo 

So coth — + 2eAI I -^^^ cos(cjt) dr, (22) 



2T 



dr 



which express Sj{uj) via A{t) or A{t) [the last equation 
is obtained using integration by parts and taking into 
account -4(0) — 0]. 

To calculate A{t) (or A), we notice that the Bloch 
equations (||)"(H) couple the dynamics of A{t) with two 
more magnitudes 



J^(r)^EnImp«2(r), 

n 

A'(T)^EnRepl'2M, 



via equations 



Ze n e 

y = ^A + ^x -ry + bimpi2, 

X ^^l:y-rX + bRcpi2, 
n 



(23) 
(24) 

(25) 
(26) 
(27) 



where b = [(/i+Za^)^/^ - (/r/2")^/2]/e. Because of the 
stationary initial conditions z(r) — pi2(t) = 0, so the 
equations are further simplified and become a closed sys- 
tem: 



A=iAI/2e)~AiH/h)y, 

y = {H/h)A+(e/h)x-ry, 

X = -{elK)y -VX. 



(28) 
(29) 
(30) 



4 



Solving these equations with the initial condition ^(0) = 
X{0) = 3^(0) = 0, one can obtain A{t) and therefore 

In the case of a symmetric qubit, e = 0, the evolution 
of X is decoupled and one can find the analytical solution 



2e 



dAjr) _ A/ 
dT 

where 
into Eq. (|2 



exp[-rr/2] 



cos fir H sin fir 

2n 



, (31) 



y/il'^ — r2/4. Substituting this expression 
we finally obtain 



eV 

Si{tj) = S'ocoth — 



f]2(A/)2r 



(32) 



It is easy to check that at zero temperature this result 
coincides with the results of Refs.t2rt2l. Notice, how- 
ever, that it does not asstmie, weakly responding detector 
(|AJ| < /) as in Refs.ll2Hl3. On the other hand, the 
derivation of Eq. (js^) assumes low-transparency QPC as 
a detector, while a much_bpjader class of linear detectors 
was considered in Refs.E9~Ej. 

With the temperature T increase the noise pedestal 
So coth(ey/2T) increases while the spectral peak around 
fi becomes lower and wider [Fig. ||(a)] because of F in- 
crease [see Eq. ([l^)]. The integral over the peak. 



[Si{lu) ~ 5*0 coth(ey/2r)] 



2tt 



(33) 



does not depend on the temperature. (As will be seen 
later, this formula remains valid for e ^ as well.) 
The peak-to-pedestal ( "signal-to- noise" ) ratio 



Sijn) - SoCoth{eV/2T) 
S'ocoth(ey/2T) 



(A/)2 



rso coth(ey/2r) 



[coth(ey/2T)]2 2(/i -I- h 



(34) 



has an upper bound equal to 4 and decreases with tem- 
perature as well as due to finite response ratio | A/|//o. In 
particular, in the case of very strong response, | A/| — 2Iq 
{I2 = 0), the upper bound for the peak-to-pedestal ratio 
is 2 instead of 4 [see Fig. ^(b)]. Overall, the effect of the 
finite response on the spectral P|eak shape is similar to 
the effect of detector nonideality.t3 

Let us emphasize that the derivation of Eq. (^) did 
not use any assumption about the magnitude of the cou- 
pling C = fi{/S.lY I SqH between the qubit and detector, 
so Eq. (^2|) remains valid even when the oscillations are 
destroyed due to strong coupling ("quantum Zeno" ef- 
fect) and replaced by a telegraph noise. The analysis of 
the fpijtCi coupling effect is completely similar to that of 
Refs£30. In particular, the quality factor of oscillations 
(disregarding the noise pedestal) is equal to 



F C 



8 {Vh 



4/n 



coth(ey/2T) ' 



(35) 
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FIG. 3. (a): Spectral density Si{u) of the detector cur- 
rent in a weakly coupled (C = 1) and weakly responding 
(AJ/Io = 0.1) regime at detector temperatures T = 0, eV/2, 
and eV for a symmetric, £ = (solid lines), and asymmetric 
qubit with e = H (dotted lines), (b): Si{uj) for a symmetric 
qubit and weakly coupled detector {C — 1) at T = for sev- 
eral response ratios AJ/Jo — 0.1, 1, 1.5, and 2. (c): the same 
as in (b) for an asymmetric qubit with e — H. 



and the transition into overdamped regime occurs at Q < 
1/2. 

For an asymmetric qubit, e 7^ 0, the analytical solution 
of Eqs. (|2^)-(|30|) is too lengthy, so it is easier to use nu- 
merical calculations and then calculate the Fourier trans- 
form ( p^ ) also numerically [Figs, ^(a) and |(c)]. In the 
next Section we will sbawi-that the result is still equivalent 
to the results of Refs.EjEj (within the common validity 
range) . 
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IV. EQUIVALENCE TO THE MASTER 
EQUATION APPROACH 

Let us remind that the master equation approach0'0 
assumes no correlation between the detector noise and 
qubit evolution and treats the oscillating part of the cur- 
rent I(t) as being proportional to the quantum operator 
z{t), so the spectral density Si{uj) should be calculated 
as 



Siiiu) 



Jdet 



4 / Kl,{T)cOs{LJT)dT, (36) 

'o 



where Sdet is the detector noise and Kz{t) = {z(t + 
T)z{t)) is the correlation function of z. In the case 
of a weakly responding detector the detector noise 
level does not depend on the qubit state and Sdet = 
S'o coth(eV^/2T). The same formula remains valid in the 
case of moderate or strong response because at high fre- 
quency 5/(00) = 2e((/+) + (/-) = 2e(/) coth(ey/2r) 
and the average curr/?nt (/) remains to be equal to /q- 
As shown in Refs.EElO, the correlation function is 
equal to the value of zM — Pii(t) — ^22(1") obtained 
from the master equation ([l0|)-(|ll|) with initial condition 
Pii(O) — 1, ^22(0) = Pi2(0) — 0. (Actually, in those pa- 
pers only the case of a weakly responding detector at zero 
temperature has been considered; however, the method 
can be easily generalized, since formally the only change 
in the master equation is a different F.) 

Let us show that Si{uj) calculated in this way coincides 
with the result obtained from Eq. ( p2|) at arbitrary qubit 
asymmetry e. For this purpose we introduce the new 
variable a = {2e/AI)A [so that we need to prove a(r) = 
Kz{t)] and from Eqs. ([2^)-(|30|) derive a new system of 
equations 



-AiH/h) y, 
~{£/h)y~ Fx, 



(37) 
(38) 
(39) 



where y = {2e/AI)y and x = {2e/AI)X. It is easy to 
see that these equations coincide with the master equa- 
tions (p^)-(pl|) for z, Impi2, and Repi2, respectively. 
Since A(0) = X{0) = 3^(0) = 0, the initial conditions for 
new variables are a(0) = 1 and a;(0) — y(0) — 0, i.e. ex- 
actly the initial conditions for Kz calculation. Therefore 
a(T) = Kz{t) and the spectral density Si{uj) calculated 
using the Bloch equations coincides with the result of the 
master equation approach. 

Hence, the analysis of Si{io) at finiia oubit asymmetry 
e is completely similar to that of Refs.E3il3. In particular, 
finite e leads to a decrease of the spectral peak around 
frequency and origination of an extra peak around zero 
frequency (Fig. ||), while the integral (|3^) does not change 
[this is a consequence of Kz{0) = 1]. An analytical ex- 
pression 



5*7(0;)= So coth 



eV £2(A/)V(4iJ^F) 

2T ^ 1 + (w n^vf^/AmTY 



(A/)V[F(l+£V2g2)] 
1+ [(c^-^!)2/[F(l + £2/a2^!2)]]■ 



(40) 



can be obtained in the limit F 51. 

Notice that for both Bloch and master equation ap- 
proaches the case of a finite detector response does not 
formally differ from the case of a weakly responding de- 
tector (only the value of the ensemble decoherence rate 
F changes). For the Bayesian approach these two cases 
are sigaificantly different, so the generalization of the 
resultalj for Si{uj) to a finite detector response (consid- 
ered in the next Section) is not trivial. 



V. GENERALIZATION OF THE BAYESIAN 
RESULTS FOR Si{uj) 

Tke Bayesian results for Si {to) have been derived in 
Ref.t3 only for the case of a weakly responding detector, 
in which the detector current I{t) can be considered as 
continuous-, and the qubit evolution is described by the 
equational^'t 



Pii 
P12 



2AI H 
-P22 = P11P22 [I{t) - /q] - 2 — Impi2, (41) 
So n 

-{pii ~ 922)-^ [I{t) ~ lo] Pi2 - 1P\2 

•JO 

. e .a , 

f 1 ^ P12 + 1 ^ (Pii - P22), 



(42) 



where 7 = F — (A/)2/45o is the qubit decoherence 
rate without ensemble averaging [for our model 7/F = 
cosh"^(ey/2r), so that 7 = at T = 0] and the statis- 
tics of l(t) can be modeled as 



m - /o = [pii(i) - P22(t)\ Ai/2 + 



(43) 



where ^(t) is the white noise with the spectral density 
S^ = So- Notice the significant difference in the mean- 
ing of pij in the Bayesian equations and in the master 
equation since Eqs. (pl|)-(p2|) describe individual qubit 
evolution without ensemble averaging. Also notice that 
the detector noise ^{t) is now significantly correlated with 
the qubit evolution pij {t) . 

To consider the case of a detector with finite response 
factor I A/ 1 /Jo, we necessarily need to take into account 
individual tunnel events in the detector because F be- 
comes comparable to lo/s- Hence, the current is not 
continuous any more and we have to use generalized 
Bayesian equationsO (which are essentially jSimilar to the 
equations of the quantum jumps formalismE3) . Then the 
qubit evolution during the time intervals between tun- 
nel events in the detector is continuous and given by the 
small-time Bloch equations for p^^ with the restored nor- 
malization: 

Pii = -P22 = -— coth(^)piip22 - 2 ^Impi2, (44) 
e 21 n 
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Pi2 = ^ coth( — ) (pii - P22) P12 

s U 
+ i-pi2 +i — {pii - P22), 



(45) 



while each tunnel event in the detector (at time t = tk) 
causes abrupt change (collapse) of the qubit state: 



/ipii(ife-O) 



Pll(ife + 0)-— n^ I r /-+ n\ 

hpil[tk - 0) + l2P22(tk - 0) 

P22(ife+0) = 1-/911 (<fc+0), 



Pl2(tfc +0) 
Pl2(tfe-0) 



Pll(^fc+0)p22(^fc+0) 
Pllitk - 0) P22(tfe - 0) 



1/2 



(46) 
(47) 

(48) 



[Actually, in Eq. ( ^ ) instead of it is better to write 
if the tunneling is in the positive direction and /~ if it 
is in the negative direction; however, the corresponding 
temperature factors cancel each other.] It is interesting 
to note that the generalized Bayesian equations do not 
contain any decoherence term even at finite temperature. 
This is because our model ai the low-transparency QPC 
describes an ideal detectom and counting tunnel events 
in both directions gives more information than measure- 
ment of only total current {I'^ — I^) assumed in Eqs. 

For the evolution simulation Eqs. (M)-(p8[) should be 
complemented by the statistics of tunnel events in the 
detector. This statistics is described by the (varying) 
rates and p~ of tunneling events in the positive and 
negative directions, respectively: 



p+(t) = (/+/e)pii(t)-K/+/e)p22(t), 
P'{t) = Ur/e)pii(i) + {I2 /e)p22{t). 



(49) 
(50) 



It is important to notice that ensemble averaging of evo- 
lution equations (H)^(48) over random moments, of tun- 
neling events described by Eqs. (K9[)-( |50|) leadsti to the 
conventional master equation (|lo|)-( ]ll])~ with the ensem- 
ble decoherence rate T given by Eq. (|12|) . 

Taxalculate Si{uj) we will use the method developed in 
Ref.Ej and write the current correlation function Ki{t) = 
{I{t + T)I{t)) = Ki{-t) at t > as 

Ki{t) = .sS{r) + (/+) e [p+ {r\+) - p- {t\+)] 

-{r)e[p+{T\-)-p-{r\-)i (51) 

where s — Si{oo)/2 determines the pedestal of Si{lo) and 
P^(t|±) is the average rate of tunneling in the positive 
{p'^) or negative {p~) direction at time t + for the 
condition that at time t a tunneling in the positive (|-t-) 
or negative (|— ) direction has occurred. 

The value of s should be chosen in a way to provide the 
correct value of Si{oo) — 6*0 coth(ey/2T) which can be 
calculated in the same manner as in the previous Section. 
For the calculation of p^(r|-|-) let us notice that as seen 
from Eq. (p6|), after the positive tunneling (at r = 0) the 
value of z = pii — P22 is equal to 



z{t + 0\+) = 



hpii — I2P22 
hpii + hp22 



(52) 



where pa are taken before the tunneling. Averaging 
z{t + 0|+) over the positive tunneling events [or, equiv- 
alently, over time with the weight factor p'^it), which 
is proportional to the denominator of Eq. (p2|)] , we get 
{z{t + 0|+)) = {hpii - I2P22)/{I), expressed via sim- 
ple averaging over time. Since (pn) = {P22) = 1/2, 
the expression can be further simplified: {z{t + 0|+)) = 
A//2(/). Similar calculation shows that (pi2(^ + 0|+)) = 

(pi2)(/l/2)l/Va>-0. 

It is sufficient to know {z{t+0\+)) and (jOi2(t+0|+)) to 
calculate p^(r I -I- ) because of the linearity of the averaged 
evolution equatio ns ([lO|)~(pT| ) in terms of z and pi2 and 
linearity of Eqs. (|49H|0). Using Eqs. (©-(P) with 



averaged initial conditions at t = 0, we can show 



p±(r|- 



AI 



2e 2(7) 



(53) 



where z(t) is calculated from Eqs. (|lO|)-(|ll|) starting 
from initial condition pu = 1, P22 = P12 = 0. Finally 
noticing that the expressions do not depend on the di- 
rection of tunneling at r = and combining the terms in 
Eq. (O) we obtain 



Kj{r>0) = {ir+^-^zir) 



(54) 



with the same z{t) as above. The constant term (7)^ 
does not contribute to the Fourier transform Si(lu) = 
2 Kj{t) cos{ujt) dr (formally it leads to a (5-function 
at cj = 0), while the second term of Eq. ( |5^ ) gives the 
same contribution as the second term of Eq. (|3^). Con- 
sequently, the calculation of Si{lu) using the generalized 
Bayesian approach leads to the same result as the master 
equation approach. 



VI. CONCLUSION 

The main result of this paper is the development of 
a method of S'/(w) calculation for a detector measuring 
quantum coherent oscillations_pf a qubit, based on the 
microscopic Bloch equations,E£l which couple the qubit 
and detector degrees of freedom. As the detector we as- 
sumed a low-transparency QPC (tunnel junction). We 
have shown that Si{uj) calculated in this wayjioraaally co- 
incides with the results obtained previouslyll3il3 by the 
Bayesian, master equation, and quantum jumps meth- 
ods, though in this paper we have considered a wider 
validity range. In particular, our formalism takes into 
account finite detector temperature T and assumes arbi- 
trary detector response |A/|//o. 

Besides that, we have generalized the Bayesian method 
of Si{uj) calculation to the case of arbitrary detector re- 
sponse and temperature. (The generalization of the mas- 
ter equation method is formally trivial.) We have proven 
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that the results of all three methods still coincide in such 
generalized case. 

The model we have considered describes essentially 
an ideal detector. The detector nonideality can be 
phenomenologically taken into account by introducing 
an extra dephasing term into the Bloch, master, and 
Bayesian equations, containing pi2. This will lead to 
an increase of the ensemble decoherence rate T and, 
therefore, to a wider and lower pesk of Si{ijj), corre- 
sponding to the qubit oscillations.E^TL^ Similar proce- 
dure can be done to take into account arqjU|bit-controlled 
change of the detector tunneling phaselljil3 [then T = 
(2/o - - (A//cos6l)2)(2e)-icoth(ey/2r) where 

9 ~ arg(M*AM)] even thougJar-auch detector is still 
ideal in the generalized sense.clil3't3 The effect of a weak 
extra coupling between the qubit and a passive finite- 
temperature environment (with temperature different 
from T)— cpi be taken into account in a way similar 
to Refs.t2l"t3. We did not consider these effects in the 
present paper because their treatment is exactly the same 
as in previous papers. 

An experimental measurement of Si{uj) and verifica- 
tion of the upper bound (< 450) for the spectral peak 
corresponding to qubit oscillations seems to be the eas- 
iest experiment related to a continuous monitoring of a 
nontrivial single qubit evolution. This makes it prefer- 
able for sooner realization in comparison with niQre |i nt er- 
esting but more difficult proposed experimentalj'c2[E3 on 
monitoring and continuous collapse of a solid-state qubit. 

The authors thank L. Fedichkin for useful discussions. 
The work was supported by NSA and ARDA under ARO 
grant DAAD19-01-1-0491. 
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